Strain-induced Fermi contour anisotropy of GaAs 2D holes 
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We report measurements of magneto-resistance commensurability peaks, induced by a square ar- 
ray of anti-dots, in GaAs (311)A two-dimensional holes as a function of applied in-plane strain. The 
data directly probe the shapes of the Fermi contours of the two spin subbands that are split thanks 
to the spin-orbit interaction and strain. The experimental results are in quantitative agreement with 
the predictions of accurate energy band calculations, and reveal that the majority spin-subband has 
a severely distorted Fermi contour whose anisotropy can be tuned with strain. 
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Manipulation of the spin-orbit coupling and the result- 
ing spin-splitting in two-dimensional (2D) carrier systems 
is of considerable current interest Such systems are 
the basis for novel spintronic devices 0, Q and also al- 
low studies of fundamental phenomena such as various 
phases that the 2D carriers acquire in mesoscopic struc- 
tures [3, 0] , and the spin Hall effect • An important 
problem is the ballistic transport of the carriers in these 
systems and, in particular, the ability to resolve and ma- 
nipulate such transport for the two spin-subbands. Mea- 
surements of the commensurability oscillations Q and, 
more recently, magnetic focusing Q, have indeed re- 
vealed spin-resolved ballistic transport in GaAs 2D hole 
systems (2DHSs). Here we report spin- resolved ballistic 
transport measurements in GaAs 2DHSs as a function of 
applied in-plane strain. Together with the results of our 
energy band calculations, the data reveal a strong distor- 
tion of the spin-subband Fermi contours with strain, es- 
pecially for the majority subband. Such distortion could 
find use in ballistic spintronic devices that rely on the 
spatial separation of carriers with different spin . 

Figure 1 captures the essence of our study. In Fig. 
1(a) we show the calculated Fermi contours of the heavy 
and light heavy-hole (HHh and HHl) spin-subbands for 
a GaAs (311)A 2DHS at different values of strain (e) 
applied along the [Oil] direction. There is a severe dis- 
tortion of the contours, especially for the HHh band, as 
a function of strain. Such distortions have been theoret- 
ically reported previously but experimental evidence has 
been only through piezo-resistance measurements which 
qualitatively agree with the calculations , [lH . In our 
study, we perforated the 2DHS with anti-dot (AD) lat- 
tices, and measured its magneto-resistance (MR) in an 
L-shaped Hall bar mesa. In the AD lattice, the carriers 
move along real-space cyclotron trajectories with orbit 
shapes that are similar to the shape of the Fermi contour 
in reciprocal space but are rotated by 90°. When these 
orbits become "commensurate" with the AD lattice pe- 
riod, the MR exhibits a peak 

The primary (smallest) such orbits are shown in Fig. 1(b) 
for the HHh and HHl holes for the current parallel to the 
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FIG. 1: (Color online) (a) Calculated Fermi contours for the 
HHh and HHl bands of GaAs (311) A 2D holes at p = 2.45 x 
10^^ cm^'^ at the indicated strain (e) values, (b) The primary 
commensurate orbits for the HHh and HHl bands for e = 
5 X 10"'' for current along [233] (see text). 



[Oil] direction. The MR peaks measured along [Oil] and 
[233] should then reveal direct information regarding the 
size and shape of the HHh and HHl contours. Here we 
provide such MR data and compare them with the peak 
positions expected from the calculated Fermi contours. 

Our sample was grown on a GaAs (311)A substrate 
by molecular beam epitaxy and contains a modulation- 
doped 2DHS confined to a GaAs/AlGaAs heterojunction. 
The Alo.35Gao.65As/GaAs interface is 100 nm below the 
surface and is separated from a 17 nm-thick Si-doped 
Alo.35Gao.65As layer by a 30 nm Alo.35Gao.65As spacer 
layer. The 2DHS has a typical low temperature mobility 
of 1 X 10^ cmVVs along [Oil] and 4.3 x 10^ cm^Vs along 
[233] at a 2D hole density, p, of 2.9 x 10" cm^^^ We fabri- 
cated an L-shaped Hall bar, with its arms along [Oil] and 
[233], via optical photo-lithography. We then deposited 
a layer of PMMA and patterned the square AD arrays 
using electron beam lithography. The AD pattern was 
etched to a depth of ~ 50 nm; this depth is sufficiently 
large to strip the dopant layer and thus deplete the holes 
in the AD regions, but small enough to avoid introduc- 
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FIG. 2: (Color online) The solid curve shows A/(= /+ - /_), 
determined from the calculated magneto-oscillations, as a 
function of strain. The circles are the experimentally mea- 
sured A/. The square is the A/ measured on a sample from 
the same wafer not mounted on a piezo. The inset shows the 
Fourier power spectra (F.P.S.) of the measured and calculated 
magneto-oscillations at e = 5.3 x 10^*. 



ing non-uniform strain in the AD lattice. In our sample 
there are two AD lattice regions, with periods a = 600 
and 750 nm, patterned on the [Oil] arm of the Hall bar, 
and two similar AD lattices, patterned on the [233] arm. 
The ratio d/a for each AD cell is ~ 1/3, where d is the 
AD diameter. The sample was thinned to ~ 200 /im 
and then glued to one side of a commercial piezo-electric 
stack actuator with the sample's [Oil] crystal direction 
aligned to the piezo's poling direction. Voltage biasing 
the piezo results in in-plane strain which is transmitted 
to the sample 1^ 13] ■ We specify strain (e) values along 
the poling direction; in the perpendicular ([233]) direc- 
tion, the strain is approximately — 0.38e [l9|. Longitu- 
dinal MR was measured, as a function of perpendicular 
magnetic field {B) aXT — 0.3 K. The data reported here 
were all taken during a single sample cool-down. 

Before presenting the commensurability MR data we 
first discuss our measurements of spin-splitting which 
we use to determine the amount of residual (built-in) 
strain in our sample. This residual strain exists because 
of the anisotropic contraction of the piezo-actuator to 
which the sample is glued, and depends on the cool- 
down of the piezo and the sample. We determine this 
strain by comparing the measured spin-splitting in our 
sample to the results of our self-consistent calculations 
of 2D hole bands, based on the 8x8 Kane Hamilto- 
nian and augmented by the strain Hamiltonian of Bir 



and Pikus [ij, |2lJ, \2f^ . They take into account the spin- 
orbit coupling due to both the structure inversion asym- 
metry of the GaAs/AlGaAs heterojunction as well as the 
bulk inversion asymmetry of the underlying zinc-blende 
structure [23|. These energy band calculations have suc- 
cessfully explained the spin-orbit induced spin-splitting 



and its de pen dence on density, electric field and strain in 
2DHSs BMlil. To make a direct comparison with the 
experimental data, we calculated the Landau fan chart as 
a function of B and determined the magneto-oscillations 
of the density of states at the Fermi energy P, [2^ . We 
then calculated the Fourier power spectrum of these os- 
cillations and obtained the frequencies /+ and /_ that 
correspond to the majority (HHh) and minority (HHl) 
spin-subbands. An example of such a spectrum is shown 
in the inset of Fig. 2 (lower trace). The difference A/ be- 
tween /-(_ and /_ is a measure of spin-splitting and can be 
directly compared to the experimentally determined A/ 
which we obtain from the Fourier transform of the mea- 
sured Shubnikov-de Haas (SdH) oscillations; the upper 
trace in Fig. 2 inset shows such a spectrum. 

In Fig. 2 we show the calculated A/ (solid curve) as a 
function of e. We emphasize that the calculations were 
done based on the sample structure and there are no fit- 
ting parameters. In this plot we have also included a 
point (open square) from SdH oscillations measured on 
a sample from the same wafer which was not mounted 
on a piezo and is therefore free of strain. This data 
point agrees well with the calculations. In Fig. 2 we also 
plot A/, determined from the SdH oscillations measured 
in the unpatterned regions of the piezo-mounted sample 
(open circles), where we have shifted all the measured 
data points horizontally by 5.7 x lO^"' in order to match 
the calculated A/. (Note that the relative strains for 
these data points are determined from our calibration of 
the piezo-actuator, and are known to better than about 
10%.) The plot of Fig. 2 implies that the sample is under 
~ 5.7 X 10^^ of residual (tensile) strain. This is some- 
what larger than the typical residual strain values in our 
experiments on other samples, but is not unreasonable. 
In the remainder of this report, we assume that this is 
indeed the residual strain in our sample. 

We now describe the results of our low-field MR mea- 
surements on the AD lattice regions of the sample. First 
we focus on the results for the a = 600 nm period re- 
gion with the current along [Oil]. The MR traces for 
this region are shown in Fig. 3(a) for different values of 
strain as indicated. We observe two strong MR peaks, 
labelled as A and B in the figure. The field positions of 
these peaks {Bp) are plotted in Fig. 3(a') for different 
values of e. Note that these strain values are based on 
Fig. 2 data, i.e., a residual strain of 5.7 x 10"'' is assumed, 
but otherwise there are no adjustable parameters. The 
data of Figs. 3(a) and (a') indicate that peak A moves 
to higher fields as the sample is further strained while 
peak B shifts to slightly smaller fields. As we discuss be- 
low, the positions of these peaks and their dependence 
on strain are in excellent agreement with what we expect 
from our calculations of our sample's Fermi contours. 

Examples of the calculated HHh and HHl spin-subband 
Fermi contours are shown in Fig. 1(a). When the cy- 
clotron orbit diameter in the direction perpendicular to 
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FIG. 3: (Color online) (a),(b) and (c): Low- field magneto-resistance traces for the indicated 2D hole density (p), anti-dot lattice 
period (a), and current direction. In each figure traces are shown for the indicated strain values from 6.9 x 10~* (top trace) 
to 4.5 X 10~* (bottom trace) in steps of 4 x 10~^. The resistance scale is for the lowest trace and the other traces are offset 
vertically for clarity. (a'),(b') and (c'): Comparison of the expected magneto-resistance peak positions (Bp) based on Fermi 
contour calculations and the positions of the experimentally observed peaks A and B. 



the current flow equals an integer multiple of the AD 
lattice period in that direction, the MR exhibits a peak 
ll[ll,[li,[ll,[ll,[l3,[i8l|. In other words, we expect MR 
peaks at fields Bp = 2hkp/eai where A;^- is the Fermi 
wavevector along the current direction and z = 1, 2, 3, . . . 
is an integer. In Fig. 1(b) we show the smallest {i — 1) 
commensurability orbits for the HHh and HHl holes that 
are expected to lead to a MR peak for the current along 
[Oil]. In Fig. 3(a') we have plotted (solid and dashed 
curves) the positions of the MR peaks for these two or- 
bits expected from the calculations. Obviously, there is 
an excellent agreement between the expected peak posi- 
tions and the experimentally observed values. This agree- 
ment is particularly remarkable considering that, except 
for assuming a residual strain of 5.7 x 10"'*, there are no 
adjustable parameters in comparing the MR peak posi- 
tions. 

We note that in our calculations the anisotropy of the 
Fermi contours emerges as follows. In hole systems the 
spin-splitting is greatly influenced by the energysepara- 
tion between the heavy- and light-hole states [l|. This 
splitting depends in turn on the direction of the in-plane 



wavevector relative to the orientation of the strain [22 1 
which results in the highly anisotropic spin-dependent 
deformation of the Fermi contours displayed in Fig. 1(a). 

In Fig. 4 we show similar data taken on the same AD 
region (a = 600 nm) as in Fig. 3, but here the MR traces 
are taken at different values of hole density while the 
strain is kept fixed at 5.7 x 10""*. Two strong MR peaks A 
and B are observed whose positions move to higher fields 
with increasing density. As seen in Fig. 4, there is again 
good quantitative agreement between the experimental 
data and the calculated positions of these MR peaks. 

Data for the AD region with period 750 nm and with 
current along [Oil] are shown in Figs. 3(b) and (b'). The 



MR traces are overall qualitatively similar to those for the 
600 nm period AD region in that they show two strong 
peaks, labelled A and B, which move in opposite direc- 
tions in field with increasing strain. The positions of 
these peaks are in good agreement with the results of the 
calculations, although peak B is observed at somewhat 
lower field values than expected. 

In Figs. 3(c) and (c') we show the data for the 750 nm 
period region with current along [233] . Again we observe 
two prominent peaks, which we label A and B, but here 
the peaks are seen at very low field values and their posi- 
tions are essentially independent of strain. Note that in 
this case the MR data probe kp along [233] for the HHh 
and HHl subbands. As seen in Fig. 1(a), in the strain 
range that our experiments probe (4.5 to 6.9 x 10~^), 
hp of the HHh and HHl bands along [233] are close in 
size and are both smaller than kp along [Oil]. This is 
refiected in the Bp curves shown in Fig. 3(c') which pre- 
dict two essentially overlapping MR peaks that occur at 
small fields and whose positions are nearly independent 
of strain. It is therefore likely that our observed peak A 
corresponds to the commensurate orbits of both HHh and 
HHl bands, and peak B is the low-field peak commonly 
observed in AD lattices with soft potentials [l4| . 

We also took data on another sample from the same 
wafer and with similar AD lattices, but not mounted on a 
piezo-actuator so that it is strain-free. Our results for this 
sample are similar to those reported by Zitzlsperger et al. 
[Tl ] where MR peaks in an un-strained 2DHS with AD 
lattices were studied, and corroborate our conclusions. 
The low-field MR traces for our un-strained sample show 
a peak at a magnetic field whose position lies between 
the two Bp theoretically expected for the HHh and HHl 
contours. Evidently we cannot resolve the two commen- 
surability peaks for the HHh and HHl bands at zero strain 
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FIG. 4: (Color online) Low-field MR traces for the 600 nm 
AD region for different densities at e = 5.7 x I0~*. Traces 
are offset, except for the bottom one. Theoretically expected 
MR peak positions are shown by filled circles on each trace. 

where the two peaks are theoretically expected to be close 
to each other. Given the relatively large width of the MR 
peaks, this is not surprising. More revealing is the obser- 
vation that the MR traces for the AD lattices patterned 
on the two Hall bar arms (along [233] and [Oil]) of our 
strain-free sample are fairly similar. This implies that the 
Fermi contours are nearly isotropic at zero strain, consis- 
tent with the calculations (see Fig. 1(a)). In contrast, the 
data in Figs. 3(b,b') and 3(c,c') are clearly very different, 
reflecting the strong anisotropy of the Fermi contours at 
large strains. 

While the overall agreement between the experimen- 
tal data and calculations shown in Figs. 3 and 4 is good 
(2^ . there are some puzzles. First, the positions of the 
lower MR peak (peak B) in most of the experimental 
data are somewhat smaller than expected from the cal- 
culations. Second, the variation (slope) of the peak A po- 
sition with strain is slightly faster than the calculations 
predict. We do not have a clear explanation for these 
discrepancies. They might be arising from the simplifica- 
tions of theoretical model that is based on the properties 
of the un-patterned 2DIIS. It is also possible that the 
cyclotron orbits are slightly distorted in the AD lattice 
regions (compared to the un-patterned regions) because 
of the non- ideal (e.g., not perfectly abrupt) potentials of 
the ADs [l3|- Third, for current along [Oil], some of the 
experimental MR traces exhibit a weak peak (labelled C) 
between peaks A and B (see Figs. 3 and 4). Note in Fig. 
1(a) that in the strain range of our measurements, kp 
along [233] for the HHh and HHl bands are quite close 
to each other. It is tempting to associate this peak C 
with a magnetic-breakdown-like phenomenon where the 
hole trajectory switches from the HHh to HHl orbit when 
these orbits come close to each other in fc-space. Sup- 



porting such a conjecture is the fact that, in the strain 
and density range where we observe peak C, the Fourier 
power spectra of the SdH oscillations also show a peak 
at the frequency /avc = (/+ + /-)/2 (see Fig. 2 inset). 
Such a peak could indicate a magnetic-breakdown [27] . 

Despite these puzzles, our results provide strong evi- 
dence that we are probing the ballistic transport in indi- 
vidual spin-subbands. They also show that we can tune 
the anisotropy of the Fermi contours, and therefore the 
cyclotron orbit trajectories, via the application of strain. 
Such tuning may find applications in ballistic spintronic 
devices, e.g., it could facilitate the separation of carriers 
in different spin subbands in magnetic focusing devices 
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